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COMPRESSED VECTOR-BASED SPECTO AT. AMAT,V5^T.g iwrpTHOP ANB,SYST gM 
5 FOR MONLINEAR RF BLOCKS 

CROSS-REFERENCE TO RELATED APPLICATION 

This application claims the benefit of the filing date of U.S. Provisional Application 
No. 60/428,432, filed November 21, 2002, entitled "Compressed Vector-Based Spectral 
10 Analysis Method and System for NonLinear RF Blocks", the disclosure of which is 
incorporated by reference herein. 

BACKGROUND OF THE INVENTION 

The present invention relates generally to simulation of signal processing systems, 
1 5 and more particularly to simulation of RF signal processing systems. 

Wireless communication circuits and systems are becoming increasingly common. 
Wireless cormnunication circuits and systems are often tested by simulation before 
fabrication to ensure correctness of operation and to help categorize wireless device 
components. Simulation of wireless communication components allows designers to 

2 0 investigate new designs and to test design changes without expending time and money in 

producing physical samples. 

Non-Unear blocks are often present in wireless communication systems. For example, 
low noise amplifiers (LNAs) are often limited by voltage sources having upper and lower 
bounds which limit the dynamic range of the LNA. Mixers are often subject to leakage from 
25 local oscillators, with self-mixing potentially occurring with the RF signal. In addition, both 
active and passive devices of all blocks may introduce non-linearities at the device level. 

Non-linearities may introduce additional spectral components in the form of undesired 
harmonics of input signals, intermodulation, and gain compression, for example. The 
undesired harmonics noiay negatively affect circuit operation. Intermodulation, particularly of 

3 0 imdesired harmonics in a multi-tone sj^em, may distort a signal being processed, as well 

increase system complexity by requiring DC offset cancellation circuitry or complex filters. 
Non-linearities also may result in reduced gain, through clipping or other effects. 

Accordingly, preferably simulation of wireless communication systems and 
compon^ts includes simulation of non-linear blocks or elements. However, simulation of 
35 non-linear elements may be difficult. Simulation of non-linearities may be extremely 
expensive from a computational perspective, with for example extraneous signal components 
spread about the continous spectrum. This computational expense is fiirther increased if the 
non-linear blocks include memory. 
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SUMMARY OF THE INVENTION 

The invention provides a simulation system and method for simulating RF systems 
5 including non-linear elements. 

One aspect of the invention provides a method of simulating radio frequency signal 
processing circuitry, comprising forming a compressed vector based equivalent of a signal; 

performing processing on the compressed vector based equivalent to simxilate radio 
frequency circuitry operation, the processing forming a processed compressed vector based 
10 equivalent of the signal; and forming an output signal using the processed compressed vector 
based equivalent of the signal. 

Another aspect of the invention provides a method of modelling circuitry, comprising 
converting first signals to compressed equivalent signals; processing the compressed 
equivalent signals to form further compressed equivalent signals; and converting the frirther 
15 compressed equivalent signals to second signals. 

Another aspect of the invention provides a system for performing RF signal 
processing modelling, the system comprising signal generator blocks forming compressed 
vector based equivalent signal representations; RF signal processing blocks processing 
compressed vector based equivalent signal representations; and conversion blocks converting 

2 0 compressed vector based equivalent signals to RF signal representations. 

These and other aspects of the invention are more ftdly understood in view of this 
disclosure. 

BRIEF DESCRIPTION OF THE FIGURES 
25 FIG. 1 illustrates an example RF signal; 

FIG. 2 illustrates the example RF signal with bands of the example RF signal 
translated to an orthogonal plane; 

FIG. 3 is a block diagram of a system in accordance with aspects of the invention; 

FIG. 4 is a flow diagram of a process in accordance with aspects of the invention; 

3 0 FIG. S is a frirther block diagram of a system in accordance with aspects of tiie 

invention; 

FIG. 6 is a block -diagram of a system simulated in accordance witii aspects of the 
invention; 

FIG. 7 is a graph of power spectral density (PSD) of an interference signal; 
3 5 FIG. 8 is a graph of a simulated output of a linear LNA in the system of FIG. 6; 

FIG. 9 is a graph of a simulated output of a linear mixer in the system of FIG. 6; 
FIG. 10 is a block diagram of a mixer simulated in accordance with aspects of the 
invention; 

FIG. 1 1 is a graph of a simulated output of a non-linear LNA of the system of FIG. 6; 
FIG. 12 is a graph of a simulated output of a non-linear mixer of the system of FIG. 6; 
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and 

FIG- 13 is a block diagram of a CVB time domain implementation of a function. 

DETAILED DESCRIPTION 

FIG. 1 illustrates an example RF signal. The example RF signal components centered 
about frequencies fb, fi, fi and f^. The signal components are non-zero in frequency bands 
having bandwidths Bq, Bi, B2, and B3, each centered about fo, fj, £2, and £3, respectively. As 
illustrated, elsewhere in the frequency spectrum the RF signal is zero. In actuaUty, noise and 
various other factors may result in non-zero signal components outside of the frequency 
bands but these signal components are often of relatively low strength. 

The RF signal may be viewed as similar to RF signals often encountered in wireless 
communication systems. For example, a wireless communication system may utilize a carrier 
frequency, represented in FIG. 1 as fi. Harmonics of the signal at the carrier frequency are 
found at frequencies fa and fa. A subharmonic, or a DC component, is located at frequency fO. 
Often, only the signal components about the carrier frequency and a few of the harmonics are 
of sufGci^t strength to be of interest. Thus, for the RF signal of FIG. 1, the signal 
components centered about fb, fi, f2, and fs may be considered frequency bands of interest. 

Accordingly, a^ects of the invention may be viewed as using only frequency bands 
of interest in simulation of a device, or of non-linear blocks of a device. In FIG. 2 the signal 
^ components in the frequency bands of interest of FIG. 1 are orthogonally transformed to 
lower frequencies. This is accomplished by defining orthogonal coordiaates in an orthogonal 
plane. The signal components in the frequency bands of interest are translated to the 
orthogonal plane. If the signal components are viewed as vectors, the signal components in 
the orthogonal plane may be viewed as compressed vectors, as the vector infonnation in the 
orthogonal plane is reduced compared to the vector information in the frequency plane. 

As an example, consider a 900 GSM 900 MHz appUcatioiL The application has a 
carrier frequency of 900 MHz. In the example, the application has a third order non-linearity, 
so that three harmonics are of interest. The total signal bandwidth is therefore 2.7 GHz. If a 
sampling frequency is 1 KHz, a vector of lengtti 2,7000,000 defines the sampled signal. If, 
however, the channel bandwidth is 6 MHz a significant portion of the vector contains 
information not of particular importance. 

Creation of vectors containing information for frequency bands of interest, for 
example in a matrix form, allows for significantly reduced, or compressed, vector size. As the 
three harmonics are of interest in the example, foiu- vectors may be formed, each defined by 
samples formed using a 1 KHz sampling frequency in the frequency bands of interest. Each 
band provides a vector of length 6,000, and the four frequency bands of interest therefore 
provide a matrix with 24,000 entries, as opposed to 2,700,000 entries otherwise provided. A 
reduced mmiber of computations may be performed in view of the reduced matrix size. 
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In aspects of the invention, sununation, multiplication, and convolution operations are 
also modijSed so as to provide operations for the translated signal components that are 
substantially eqmvalent to those operations for untranslated signals. In aspects of the 
invention, these operations are defined as follows. 

A discrete signal x[n] is a Piecewise Non-Zero (PWNZ) signal if it can be described 
as follows: 

E.W 
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A discrete signal x[n] is a Equally Piecewise Non-Zero (EPWNZ) signal if it is PWNZ with 
the following properties: 
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where Lp is the piece length 
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Using the compressed vector terminology of above, a compressed vector-based (CVB) 
equivalent of an EPWNZ signal is: 



or Xyp = 



Where x, e /e'"^' and Xype i?''"^' 



(3) 



Summation of CVB equivalent signals may be considered as: 



35 



Given Zvp[n] as the summation of Xvp[n] and yvp[n] (i.e. Zvp[n] = Xvp[n] +vp yvp[n] where "+vp" 
is a symbol of sunmation of two CVB equivalent signals), then 
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Where x^.y^^andz^ e JJ**"^' and .yyp.andzyp e J?^'^'. Also Vis the component wise 



summationof twolxLp vectors. 



Subtraction of CVB equivalent signals may be considered as: 
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Given Zvp[n] as the subtraction of Xvp[n] and yvp[n] (i.e. Zvp[n] == Xvp[n] yvp[n] where is a 
symbol of subtraction of two CVB eqxiivalent signals), then 
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WhereXf , , andzj e and Xyp , yvp, andzyp € R^"^"" • Also"-"is thecomponentovise 
subtractioEiof twolxLp vectors. 
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Multiplication of CVB equivalent signals maybe considered as: 
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Given Zvp[n] as the multiplication of Xvp[n] and yvp[n] (i.e. Zvp[n] = Xvp[n] Xyp yvp[n] where 
"Xyp" is a symbol of multiplication of two CVB equivalent signals), then 



Zyp[n] = Xyp[n]xyyp[n] 

VP 



















35 









(6) 



J^erexpj/^andz. ^R^"^^ and Xyp ,y^,2cndzy^ eiJ^*^'. Also Vis the componentwise 
multiplicaionoi twolxLp vectors. 



Division of CVB equivsdent signals may be considered as: 
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Given Zvp[n] as the division of Xvp[n] and yvp[n] (i.e. Zvp[n] = Xvp[n] -ryp yvp[n] where "-^vp" is 
a symbol of division of two CVB equivalent signals), then 
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Wherexj.y^.andz^ ei?"^ and , yyp, aadzyp €i?^'^''.Also"V4sthecomponen1wise 
divisionof twolxLp vectors. 

Scalar multiplication of a CVB eqxiivalent signal may be considered as: 

Given Zvp[n] as the scalar multiplication of Xvp[n] with the scalar d (i.e. Zvp[n] = d.vpXvp[n] 
where ".vp" is a symbol of a scalar multiplication of a CVB equivalent signal with a scalar), 
then 
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Wherex^.BXidz^^ g B}""^^ and Xyp ,andzvp e i?^"^^. Also"m:/4s the componentwise 
scalar multiplic^on of al x Lp vector and a scalar. 

Convolution of CVB equivalent signals may be considered as: 

30 

Given 2vp[n] as the convolution of Xvp[n] and yvp[n] (i.e. Zypln]^ Xyplnlf^yypin] where 

VP 

"(g)" then 

VP 
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i=2M-l 
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In which xr and ypj are flipped version of Xj and yi reqjectively (i.e. 
Xfi(^0 = Xf(length(X()-h)). M is the maximum number of sub-pieces taken into account, 

and 0 is as follows: 

X, (^yj = Jx,(F)^/F -f)dF (10) 

—00 

In some aspects of the invention, x,<^yj is substituted with a regular convolution if 

yj(F-f)=yjif~F). 

FIG. 3 is a block diagram of a simulation system vising frequency bands of interest, in 
some aspects in terms of compressed vectors as discussed above. In the block diagram of 
FIG. 3, an input signal 301 (or signals) is provided to a converter 303. The converter converts 
the input signal to a CVB signal. The conversion in some embodiments makes use of 
appropriate ^roximation functions such as a sine function or otiier functions. The converter 
provides the CVB signal to a CVB processing block 305, which perfomis various operations 
as described on the input signal and signals generated internal to the CVB processing block. 
The signals generated internal to the CVB processing block may be formed by internal signal 
generation blocks, by processing the input signal, by various combinations of the two, or 
signals ultimately generated from a combination of some or all of these signals. The CVB 
processing block provides resulting signals to a reconversion block 307. The reconversion 
block converts the signals back to the form of the input signals, either time domain signals or 
frequency domain signals as the case may be. If appropriate. Fast Fourier Transform (FFT) 
and inverse FFT (iFFT) processing may be used. The reconversion block provides an output 
signal 309 (or signals). 

In various aspects of the invention the CVB signals are m terms of dB, while in other 
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aspects of the invention other scales are used. 

FIG. 4 is a flow diagram of a process of simulating a component or system using 

5 frequency bands of interest. In Block 401 an input signal (or signals) is converted to a CVB 
format. In block 403 the system is simulated using CVB operations, for example as described 
above, and CVB format resxilts are produced. In block 405 the CVB format results are 
reconverted to the form of the input signal and provided as an output. 

FIG. 5 is a block diagram of a CVB processing block. The CVB processing block 

10 receives an input signal 501 in CVB format. The input signal is provided to a first simulation 
block 503. The first simulation block performs operations as described above and below. The 
CVB processing block also includes a first signal generator block 505 and a second signal 
generator block 507. In some embodiments, the signal generator blocks are white Gaussian 
noise generators, adjustable single tone generators, adjustable two tone graerators, adjustable 

15 mult-tone generators, GSM900 standard power spectral density (PSD) mask generators, 
GSM900 adjustable PSD mask generators, GSM standard Interferer generators, GSM 
adjustable Interferer generators, or adjustable white Gaussian noise generators or oscillators. 
The signal generation blocks may be developed using the operations described above and 
below. 

20 In the CVB processing block of FIG. 5 the first signal generator block provides a 

signal to a second simulation block 509. The second signal generator provides a signal to the 
first simulation block and the second simulation block. The first simulation block performs 
operations as described above and below using its input signals. Similarly, the second 
simulation block performs operations as described above and below using its input signals. 

25 The first simulation block and the second simulation block each provide a signal, formed by 
the performing of operations, to a third simulation block 511. The third simulation block 
performs operations as described above and below using its input signals. The third 
simulation block provides an output signal, formed by the performing of operations. 

In various embodiments the simulation blocks may be considered as combinations of 

30 the operations described above. For example, the simulation blocks may be Linear Time 
Invariant without memory (LTI) blocks. Linear Time Invariant with memory (LTIM) blocks, 
Non-Linear Time Invariant without memory (NLTI) blocks, and Non-Linear Time Livariant 
with memory (NLTIM) blocks. These blocks may in turn be used to form low noise amplifier 
blocks, mixer blocks, and other blocks. 

35 As a notation convenience, an input to a block may be considered to be x(t) and an 

output of a block may be considered y(t). Also as a notation convenience X(f) and Y(f) may 
be considered Foiirier tranforms of x(t) and y(t), respectively. 

For an LTI operation y(t)=dix(t) and Y(f)= diX(f), where di is a constant scalar. The 
CVB equivalent representation of the LTI operation, the LTI block, in the fi-equency domain 
isYvp(f)=di.vpXvp(f). 
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For an LITM operation y(t) = jhiT)x(t-T)dT and Y(f)=H(f)X(0, where H(f) is the 

—CO 

5 Fourier transform of the impiilse response h(r) . 

The CVB equivalent representation of the LTIM operation, the LTIM block, in the frequency 
domain is 7^ (/) = i/^ (/) x ^vp (/) , where Hvp(f) is the CVB equivalent of H(Q. 

VP 

For an NLTI operation, 

10 

y(t) = ao + a, jp(0 + a jjc^ (0 + a^x^ (0 + ... 

Where 0^1,0:2,... are coefficients representing non-linearity of a device and * indicates 
convolution of two signals. The coefficients may be extracted, for example, by performing 
curve fitting on a curve relating inputs and outputs of the device. The CVB equivalent 
representation of the NLTI operation, the NLTI block, in the frequency domain is 



VP VP VP VP VP VP 

20 

In some embodiments NLTIM blocks are modeled as combinations of NLTI sub- 
blocks and LTIM sub-blocks. In other embodiments a Volterra series method is used. In a 
Volterra series an n* order non-linear operator is described as 

00 00 00 

2 5 y(.0 = jfh (t, )x(t - r, + ... + J... \h„ (r, ,..r„ )x(t - r, ). . jc(f - r„ )dT, ...dT„ 

-«0 —CO —00 

in which x(t) is the input, y(t) is the output, and (ri„..T J = 0 for any < 0 y = l,...,n ), 
Further information may be found in Schetzen, M., "The Voleterra and Wiraer Theories of 

2 Q Nonlinear Systems" (1 980), the disclosure of which is incorporated by reference. 

Simulation may be computationally expensive when convolution operations are 
performed in the frequency domain. Accordingly, in some embodiments of the invention 
CVB equivalent signals in the time domain are used, allowing for, for example, 
multiplication in the time domain instead of convolution in the frequency domain. A real time 

25 domain signal s(t) has a Fourier transform S(f). Svb(f) is a CVB equivalent representation in 
the frequency domain. The CVB equivalent representation in the time domam is Svb(t), with 
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and 



10 



SyBt (0 = [5(0.exp(-72;z^c0] * [(■;^).sin(;zSO] 

Ttt 



which utilizes an. i^roximation function, and where fc is the the carrier frequency and B is 
ttie bandwidth of interest For the fiequency domain 
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with 



25 



CVB equivalent systems are formed using a combination of frequency shift blocks 
and samplings. For example f(x)=x^. 
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fvBl i^VB ) 



may be implemented as shown in FIG, 13, with a sampling frequency greater than the signal 
bandwidth. 

FIG. 6 is a block diagram of a simplified RF receiver front-end. The front-end 
includes an LNA 601 receiving an input signal. The LNA ampUfies the input signal to 
provide a first input to a mixer 603. The mixer receives as a second input the output of a local 
oscillator (LO) 605. The mixer mixes the two signals to form an output signal. 

As an example, the input signal is a GSM900 standard Interference signal. FIG. 7 is a 
graph of the input signal as provided by a CVB signal generation block. FIG. 8 is a graph of 
an output of the LNA when the LNA is provided the input signal of FIG. 7 and the LNA is 
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modeled as an LTI block. Similarly, FIG. 9 is a graph of an output of the mixer of FIG. 6 
when the mixer is modeled as a linear mixer. 
5 The system of FIG. 6 may also be modeled as a non-linear system. For example, the 

LNA may be modeled as NLTI block. In one embodiment, the LNA is modeled as an NLTI 
with 

10 With the non-linear coefficients derived from curve fitting, with curves generated for 
example using tools from Cadence corporation. In the example described the coefficients are 

a© = 0 ai = 1 0.27 =-4.5 . 

Similarly, the mixer may be modeled as a non-linear block, hi the example described, the 
15 mixer is modeled m accordance with the block diagram of FIG. 10. In FIG. 10, the mixer is 
an NLTI block 1001. The NLTI block receives an input fix)m an LO 1003. In addition, with 
RF input to the mixer, LO to RF leakage is modeled using a gain block 1005, with LO to RF 
isolation being 30dB. RF to IF leakage and RF to LO leakage are assimied negligible. The 
mixer is modeled by 

20 

The coefficients are derived using curve fitting, Avith for example curves generated using 
tools from Cadence Corporation. In the described example, the coefficients are 

25 

ao = 0 flfj = 1 = 0.32 ^3 = -3.21 . 

Simulation results are shown in the graphs of FIGs. 1 1 and 12, with the graph of FIG. 
11 shoAving the LNA output and the gK^h of FIG. 12 showing the mixer output. The 
30 simulation time, using an Intel Corporation Pentium 4 equiped personal computer wifli a 
fi^uency resolution of 5KHz and a 6GHz signal bandwith, was less than one minute. 

Accordingly, the invention provides a modeling system and method using compressed 
vectors. Although the invention has been described in certain embodiments, it should be 
recognized that the invention encompasses the claims supported by this disclosure and the 
3 5 their equivalents. 
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